Transition Probabilities for 
Telephone Traffic* 

By V. E. BENES 

(Manuscript received April 21, 1960) 

-4 stochastic model for the occupancy N(t) of a telephone trunk group is 
specified by the conditions that arriving calls form a renewal process, that 
holding times have a negative exponential distribution, and that lost calls are 
cleared. The transition probabilities of N(t) are determined, and their limits 
are studied. These transition probabilities have practical value in making 
theoretical estimates of sampling error in traffic measurements, and in the 
study of overflow traffic. 

I. INTRODUCTION 

We shall study a stochastic process {N(t), t ^ 0), which is a mathe- 
matical model for the occupancy of N service facilities, with no provisions 
for delays. For example, N(t) can be interpreted as the number of (fully 
accessible) telephone channels (trunks) out of a group of N such in use 
at time /, with lost calls cleared. Also, we can think of N(t) as the num- 
ber of items on order at time t in an idealized inventory situation in 
which at most N items can be on order at one time (see Arrow, Karlin 
and Scarf ). Throughout the paper we use terminology appropriate to an 
application to telephone trunking. The process N(t) is determined by 
the following assumptions: 

i. Holding times of trunks are independent, each with the same nega- 
tive exponential distribution function, of mean y~ l , y being the "hang-up 
rate." 

ii. Times between successive attempts to place a call (interarrival 
times) are independent; each has the distribution function A( • ), where 
A(-) is arbitrary except for the condition A(0) = 0. This assumption 
covers Poisson arrivals as a special case. The mean of A(-), when it 
exists, is denoted by ^ . 

* This work was completed while the author was visiting lecturer at Dartmouth 
College, 1959-60. 
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iii. There are JV trunks, N being finite. 

iv. Calls that find all N trunks busy are lost, and are cleared from the 
system without effect on the flow of arrivals (no retrials). These or 
similar assumptions appear in Palm, 2 and in Pollaczek; 8 ' ' ,6 certain prop- 
erties of N(t) itself have been studied by Takacs, 7,8,9 Cohen 10 and BeneS." 

II. SUMMARY 

The random process of interest is N(t), which is interpreted as the 
number of trunks in use, or the number of calls in progress, at time t; 
N(t) is a random step function fluctuating in unit steps from to N. 
For the most part, we restrict attention to that version of N(t), written 
N(t — 0), that is continuous from the left. 

The present paper is chiefly theoretical in character. It provides (a) 
formulas for the Laplace transforms of the transition probabilities of the 
stochastic process N(t - 0), and (b) a statistical description of the calls 
that overflow a trunk group of the kind described in Section I. The for- 
mulas will be exemplified and used in a second paper, 12 where specific 
applications to switch counting and traffic averaging are described. 

We begin Section III with a general account telling what transition 
probabilities are and why they are useful and interesting in traffic 
theory. The primary result, Theorem 1, can then be stated; it completely 
characterizes the transition probabilities 

Pr {N(t - 0) = n|AT(0+) = m] 

as functions of t by determining their Laplace transforms, under the 
restriction that A ( • ) has a probability density. Section III ends with 
a computation of some important transition probabilities for Poisson 
arrivals; practical consequences of these results will be developed in the 
second paper. 

We prove Theorem 1 in the Appendix A. If y(t) is the time elapsed 
since the last call arrival prior to t, the process [N(t - 0),y(l)\ is Mar- 
kov, and we calculate its distributions from the usual Kolmogorov equa- 
tions. The stationary distribution of this Markov process is determined 
in Appendix B. 

In Appendix C the process N(t - 0) is studied directly in terms of 
renewal theory and regenerative processes, using results of Smith. No 
assumptions of absolute continuity arc made. This procedure leads to an 
extension of Theorem 1, and other results outlined in the next para- 
graphs. (Details are omitted.) 

Let R n be this event: a call arrives and finds n trunks in use. Each 
occurrence of R n , where g n ^ N, is a regeneration point oiN(t -0), 
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in the sense that the history of N(t — 0) prior to the given occurrence 
of R n is statistically irrelevant to the development of N(t — 0) after 
the occurrence. Let x m ,„ be the time elapsing from an occurrence of R, n 
to the next occurrence of R„ ■ We prove x m ,„ < co with probability one, 
and, if 



Mi 



= / xdA(x) < =o , 



then E{x m , n \ < oc. 

The underlying probability functions that we calculate in Appendix 
Care, for ^ n ^ N: 

00 

Q n (t) — 22 Pr U>'th call arrives before t and finds n trunks busy) 

(1) 
= #{ number of occurrences of R n in [0,0} • 

From this interpretation it is apparent that the (?„(•) are unbounded 
monotone functions; one may expect them to be ultimately linear. The 
transition probabilities of N(t — 0) can be represented in terms of the 
functions Q„(-) and the transition probabilities of the simple death 
process with death rate 7 per head of population, if the Q n ( ■ ) are eval- 
uated for appropriate initial conditions. This is done in Appendix C. 
With this representation we investigate the existence of 

limPr{AT(J - 0) = n). 

(-00 

From Theorem 4 and the solutions for the Laplace-Stieltjes transforms 
of the Q„{ ■ ), this limit, when it exists, can be evaluated explicitly, using 
the relation 

Pn 

where p„ is the equilibrium probability that an arriving cell will find n 
trunks in use. (For p n see Refs. 7 and 11.) 

III. TRANSITION PROBABILITIES OF N(t) 

The transition probabilities of a stochastic process x t tell how likely it 
is that the random function .r ( ., take on a value z at a time /, if it is 
known that it took on the value ij at time s. Such a transition probability 
is written 

Pr{.r, - z\x s = y\, (2) 

the vertical bar being read and interpreted as "given that" or "if." 



(3) 
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In other words, (2) expresses the relevance of the information that the 
event {.r„ = y\ has occurred to the likelihood that the event {x t = z\ 
will occur. In still other words, (2) expresses the dependence of the 
event [x t = z) on {x 8 = y\. 

The chief practical use of transition probabilities for models of tele- 
phone traffic is in computation of covariance functions; these, in turn, are 
used to compute theoretical estimates of sampling error in actual traffic 
measurements, such as time averages and switch counts. To sec how 
this happens in a particular case, we consider the use of the continuous 
time average 

M(T) = % f N(t)dt 

as an estimate of the carried load. The variance of M is 

E{M 2 \ - E-[M\ = 

T* C [ T [E[N(t)N(s)\ - E{N(l)\E{N(s)}]dsdt. 

Jo Jo 

The integrand is just the covariance R(t,s) between N(t) and N(s); 
if N(t) is stationary in the wide sense, so that R(t,s) = R(t — s), then 
(3) reduces (by partial integrations) to 

C T 
Var {M} = 2T' 2 (T - t)R(t) dt. (4) 

Jo 

The covariance R(t) can be written in terms of the transition proba- 
bilities of N( ■ ) as 

N N /N \2 

fl(i) = IZ mnp m Pr \N(t) = n \ N(0) = m) - ( £ mp m ) , (5) 

m=0 n=0 \m-0 / 

where [p m ] is the stationary distribution of N( ■ ). Formulas (4) and (5) 
then indicate how the transition probabilities can be used to find the 
variance of M . 

Our basic result concerning transition probabilities is most easily ex- 
plained and understood after some of the notions used in stating it are 
discussed. The first few are merely abbreviations; we let 

A*(s) = f° e-'dA(t) = a (s), 
Jo 

a n (s) = A*(s 4- ny) , 

X = 1, 

v 1 - a n (s) v W 1 - aj (s) 

A n -7—7 -A n _i — ^J_ / \ • 

a n (s) /-i aj{s) 
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These Laplace transforms enter because we shall be characterizing the 
Laplace transforms of the transition probabilities in terms of the hang-up 
rate y and the transform A*(s) of the interarrival probability density. 
In the summary we have denoted by R n the event: a call arrives and 
finds n trunks in use. We let q n (t,0) be the "density" of R n at time t, 
that is, the rate at which R n is occurring at t, and we let 

WO - t({)qj(t,0) (6) 

be the associated binomial moment. From (1), it can be seen that 
dQn/dt = q„(t,0), when the former exists. The b n ( • ) and the q H ( -,0) are 
also related by the inversion formula 

q n (t,0) = Z (-D'(* J" J W). 

j=o \ n / 

More generally, we use q n (t,u) as a density function in the variable u 
with the heuristic meaning 

q„(t,u)du = Pr [N(t - 0) = n and u < y(l) ^ u + da). 

We can now state 

Theorem 1: The transition probabilities of N(t — 0) may be deter- 
mined from the generating function formula 

E{x N(t - 0) \ = fj2q H (t - 2/,0)[P,(s)] B+1 -*»[l - A(y)]dy 

(7) 

where 

I\(x) = 1 + (x - l)e- yu , 

and the Laplace transforms of the binomial moments &«(•) are given 
by 

b,*(s) = (X,,)-' (bo* - t [(j H. x ) *** ~ ' 



an(s)_ 



AVi , 



h*(s) - 



1 - A*(s)' 



fco* , y^ Xj-\kj*(&) 



»=o \n / 
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where 

„ f A /j\ ,_ s - nyt a(l + u) du 
k n * = Laplace transform of j L ( n ) </;(0,w) e ^ _ j^t . 

The k„* introduce dependence on the boundary conditions at t = 
expressed by the functions q n (0,u). The Kronecker symbol 8 Hn in (7) 
indicates that a call is lost if it finds all trunks busy. 

To show how Theorem 1 can be used we shall compute the Laplace 
transforms of 

Pr \N(t) = N | N(0) = m\, m = 0, 1, • • • , N, 

in the important special case of Poisson arrivals at rate a, for which a 
great simplification of the formulas occurs. In this case, 

^) = \o , /<o. 

Also, we set 7 = hang-up rate = 1, which amounts to measuring time 
in units of mean holding time; then 

o,(«)- _ L °_ L . (8) 

a + s + n 

Our choice of the transition probability to the "all trunks busy" con- 
dition \N(t) = N\ as an example is not arbitrary; it turns out that, in 
many cases, including Poisson arrivals, the mean of N(t) and the co- 
variance depend only on the transition probability to the "boundary" 
condition \N(l) = N] . A similar situation occurs in the theory of queues 
with one server: the mean delay can be written as an integral of the 
probability of being on the "boundary," i.e., the chance that the server 
is idle. 14 

Since arrivals are Poisson, the y( ■ ) process is in fact superfluous, and 
we may assume N(0) = m,y(0) = 0, so that 



*..(,) = />—«(:)«*-(:) 



n ^ m, 



a + s + n ' (9) 

= n > m. 

In formula (7) (Theorem 1) set x = 1 + w, and take Laplace transforms 
with respect to t; the coefficient of w N is 

f <T" Pr \N(t) = N I N(0) = m\ dt 



= f e - 8 <-»'-' d t[q N * + q N -i* + 8 Nm ], 
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where q H *(s) is the transform of q„(t,0). Now, from (6), (9) and (24) 
we find 

* . * 9-v* — 5 lVm a iV (s) 
g.v -t- g.v-i — -t-v ; 

hence 

f e-'Pr {#(/) = N\N(0) = m] dt = ^- . (10) 

^o a 

This result can also be obtained heuristically as follows: 

a = total density of arrivals at t 

= q M (tfl) + a[l - Pr \N(t - 0) = N\)\ 

taking Laplace transforms, we get, (10). 
From Theorem 1 and (10), we find 

a 

/m\ /m\ 1 - aais) ..... ( m \ [1 - ao(s)] •••[!- a„,-i(s) ] 
. i\0/ \1 / a a (s) ' \m/ ao(s) • • • a m _i(s) 



I — o (s) 



EfTx. 



Oo(s) 

But, for our example, (8) implies 

1 — a„(s) _ s + n 
a„(s) a 

hence, defining (after J. Riordan in the Appendix to Wilkinson 10 ) the 
"sigma" functions a k (m) by 

»' w /i i - i \ m—i _l 

with m (but not k) an integer, we can show that a~ l q N * reduces to 

[~ c-'Pr\N(t) = N\N(0) = m\dt = C^'AT? • (11 > 

This and similar results for Poisson arrivals have been found by S. O. 
Rice in unpublished work. 

t The "sigma" functions are related (o the Poisson-Charlier polynomials p„(x) = 
a"' 2 (n!)» J2 f-l)"~' ("•) j\arl (j) by a k (m) = (-o*)"(m!)-*p M (-ifc). See Szegd. 16 
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Since the event {N(t) = N] (the "all trunks busy" condition) is of 
primary interest, the transition probability 

p„(f) - Pr [N(t) = N | JV(0) = N\ 

has been used (e.g. by Kosten 17 ) as a "recovery" or "relaxation" function 
that is characteristic of the dynamic behavior of the system, especially 
of its approach to equilibrium from the "all trunks busy" condition. 
Such a function has been computed from (11) and plotted as Fig. 1, 
for a (heavy) load of 10 erlangs offered to 5 trunks, giving a loss proba- 
bility of 0.563. 



IV. OVERFLOW TRAFFIC 

In the design and engineering of trunking plans in telephony, it is 
common practice to offer the calls lost by one trunk group to a second 
or overflow group. It has been discovered that the right choices of group 
size and the pooling of overflow traffic can lead to efficient trunking 
arrangements, called graded multiples. For this reason, some theoretical 
work, as well as much empirical study, has been devoted to the statisti- 
cal behavior of overflowing calls. The principal references are to Brock- 
meyer, 18 Cohen, 10 Kosten, 19 Palm, 2 Takacs, 7,8,9 and Wilkinson. 15 

In accordance with current usage in mathematical literature, let us 
refer to a sequence of mutually independent, identically distributed, 
positive random variables as a renewal process. The interarrival times 
that we have assumed in the model describing the trunk group then form 
a renewal process. It has been shown by Palm 2 that, if calls arriving in 

i.o 



o.s 
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Fig. 1 — "Recovery function" Pr \N(t) = N | iV(0) = N] for N = 5 trunks and 
a = 10 erlangs (heavy traffic). 
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a renewal process are served by a finite group of trunks, with exponential 
holding times and lost calls cleared, then the overflowing calls can also 
be described by a renewal process. That is, the time intervals between 
successive overflowing calls are mutually independent and identically 
distributed. Palm also showed how the distribution function of these 
interoverflow times can be calculated from the interarrival distribution, 
the hang-up rate and the group size. 

We can deduce Palm's results in a simple way from our basic theorem 
and give a general formula for the Laplace-Stieltjes transform of the 
interoverflow distribution. Let 0\(t) be the average number of overflows 
occurring in the closed interval [0,/], assuming that an overflow occurred 
at time 0. Thus N (t) is the particular form of Q N (l) that arises when 
Mi = and JV(0— ) = N. We use G(u) to denote the distribution func- 
tion of the interoverflow limes. Since these times are independent, it 
can be seen that 

N (t) = U(t) + f OAt - u)dG(n), t ^ 0, (12) 

Jo 

where U(t) is 1 for t ^ 0, and otherwise. If N *(s) and G*(s) are the 
respective Laplace-Stieltjes transforms of Oat and G, then (12) implies 

G*( S ) - °' v * (s) - 1 

( ' W " N *(s) ' 

which determines G(u) uniquely if On* is known. 

Since, as noted, On* is the particular case of Q N * arising when Wj = 
and iV(0— ) = JV, a formula for it can be found from (32). In the particu- 
lar case being considered 



K* = 



©• 



and so On* is given by 

1 



1 - A*(s) 

(N\ 1 - a n (s) (N\ [l-fl ( S )]---[l-fl,v- 1 (s)] ( 13) 

^Vl/ a,(s) "*" ' ' ^\NJ ai (s) ■■■ cin(s) 

- ai(s)] •••[!- a„(s)J 



UT 



ai(s) • • • a n (s) 

If n\ = / x dA(x) < co , the mean time between overflows is 

Jo 



Mi 
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where p N is the equilibrium probability of loss (studied in Refs. 7 and 
11). 
For N = 1, (13) gives 

1 - A*(s) + A*(y + s) (u) 

0l 1 - A*(s) 

G* = A *^ + s) . (15) 

1 - A*(s) + A*(y + s) 

Since A*(y + s) is the Laplace-Stieltjes transform of 
Ll (t) = f e~ yu dA(u) , 

■ 

(15) can be inverted to give 

G(t) = X>„(0 = {5(4)1(0, (16) 

where 

f>i — Li t 

<Pn+l = <Pn * (A — Li) 

and "ir" denotes Stieltjes convolution. 

Formula (13) agrees with the recurrence relation given by Palm 2 for 
the overflow distribution from N trunks. The "one-trunk" case of (14) 
through (16) is important theoretically because all other cases can be 
obtained from it by iteration. Formula ( 16) defines a mapping G = %(A ) 
and the interoverflow distribution for N trunks can be written as 5^(4 ), 
the Nth iterate. 

For one trunk, the first two moments of the interoverflow time u are 

EM -*, 

2 2 

2, M , *f*i 



E{u) = =■ + 

(Li CL\ 



1 - - f° le~ yt dA(t)~\, 
Mi • o J 

where m = I u dA. In particular, the ratio of second to first moment 



is 



U^i = « + 2«r 1 _L/" fe -. dA (ol, 

£i\U\ Mi Oi _ Mi -n J 



TRANSITION PROBABILITIES FOR TELEPHONE TRAFFIC 1307 

so that the mapping $ always increases this ratio, by an amount pro- 
portional to E[u\. 

APPENDIX A 

Approach Using a Markov Process 

Let N{t) be the number of trunks in use at time t. To study the dis- 
tribution of N{t) we introduce the two-dimensional process \N(t),y{t)\, 
where y(t) is the length of the time interval from / back to the last 
arrival epoch prior to /. We assume that -4( • ) is absolutely continuous, 
with a continuous density a( • ). 

The reason for using the two-dimensional variate is that, unless arriv- 
als are Poisson, the N(t) process by itself is not Markov. To avail our- 
selves of the functional equations satisfied by the distributions of Markov- 
processes, we include y(t) in the "state of the system." This inclusion 
does result in a Markov process. The device of "Markovization" by the 
inclusion of variables has been suggested and developed by Cox, 20 and 
also has been used by Takacs.'' 8,9 

It is natural physically to think of the random functions N(t) and 
y{t) as being continuous from the right. However, we shall assume only 
that y(t) is always defined to be equal to y(t + 0), and shall study the 
two processes {N(t + 0),y(t + 0)| and {N(t - 0),y(l + 0)| jointly. 

That N(t — 0) and N(t + 0) are not the same process is clear: 
N(t- 0) = N and y(t) = imply N(t + 0) = A^; but, if N(t + 0) = N, 
y(t) = 0, then N(t — 0) = N or .V — 1 according as the call that just 
arrived is lost or accommodated. The analysis of N(t — 0) and ,V(/ + 0) 
shall be carried out in terms of two sets of probability density functions, 
p„(t,y) and q„{t,y), where 

p„(t,y)dy = Pr {N(t + 0) = n and y < y(t) ^ y + dy\, 

qn(t,y)dy - Pr [N(t - 0) = n and y < y(i) ^ y + dy\. 

Lemma: p„(l,y) = (],,((,!/) for almost all y. 

Proof: Let P be a basic probability measure determined by our as- 
sumptions (i) through (iv) of Section I; P is defined for sets of elements 
o) in a space £2. We assume further that N(t,co) is separable, so that 

S ( = n [e>:N(t - u) = N(t + 0) = N(t + u)\ 
is a measurable set. 
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Now if 2/(0 > 8 > e, then y(t - e) = y(t) - e, almost surely, arid 

Pr (ft \ Nit + 0), y(t) > 8) ^ <T 2 ^ l ~ ^^ 

independently of N(t + 0), almost everywhere, so that 

Pr {& 1 2/(0 > 8} ^ e ___^— . 

The sets *S £ are monotone nondecreasing, so S = lim S« as e — > is 
measurable, and 

Pr {So | 7/(0 > 5} = 1, almost everywhere [P], (17) 

and So is the w-set on which N(t) is constant in some interval (t — u, 
t-\-u). The lemma follows from (17). 

It remains to establish the relationship between p„(t,y) and q n (t,y) 
when y = 0. From our previous remarks about N(t — 0) and N(t + 0) 
it can be seen that 

p N (l,0) = q N (tfi) + q»-i(t,Q), 

Pnitfl) = q n -i(t,0), 1 <n ^N - 1, 

Mtfi) = o. 

To formulate the Kolmogorov equations for p„(t,-) and q n (t,-), we 
need the function X( • ) denned by 

This is the probability density that an interarrival time will end in the 
next dy, given that it has lasted a time y to date. The functions q n (t, •), 
where £ n 2g N, (with q N+i = 0), satisfy the difference-differential 
system 

|~£ + i. + 77l + X(y)l g. = 7(n + l)frn, (18) 

|_d* d?/ J 

and the behavior of the densities q n (t, • ) for y = is determined by the 
additional condition 

Qnitfi) = P q n (t,y)\(y)dy. (19) 

Jo 
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We introduce the generating function 

N 

M*,t,y) = Zx n qnd,y), 

n-0 

and from (18) obtain 

K+^+^-^s+H^ ' (2o) 

whose general solution is 

f(x,t,y) - K[t - y, r w (* - 1)|[1 - A(y)\. 

Before continuing, we note that the functions p„(t,y) also satisfy the 
system (18), but that, the analog of (19) is 

Pn+i-tyjtfl) = I p n (t,y)\(y) dy , (21) 

Jo 

where the Kronccker 8 symbol is used to indicate that an arriving call 
is lost if it finds all N trunks busy. The generating function <p(x,t,y) of 
the p„(t,y) is also a solution of (20). 

The function \p(x,t, • ) is ^-continuous for y > 0, so, from the lemma 
proved previously, we conclude that \f/(x,t,y) = <p{x,t,y) almost every- 
where in y, and that 

lim $(x,t,y) = <p(x,1,Q). 

Because of the "lost calls cleared" assumption, we must have 

^(.r,/,0+) = x+(x,t,0) - x*(x - l)qs(t,0) 
= <p(x,t,0), 

so that ip is discontinuous in y at y = 0. 

Let P u = P y (x) abbreviate 1 + (x — 1)<T™. It can be verified that, 
the function K(-,-) in the solution of (20) is given by 

K(u,z) = (1 + *)*(1 + z,u,0) - z(\ + z) N q N {u,0), I ^ y, 

*(1 + ze^Q-u) 
~ i-A(-u) ' "' 

for the solution \J/(x,t,y). From this we find that, for t ^ y, 

f(x,t,y) = 

PMP u ,t ~ Z/,0)[l - A(y)} - (P u - l)P„V(< - ?/,0)[l - A(y)], 
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while, for t < y, 

+(x,t,y) = f(P<,0,V - t) 1 1 _' A fy y l t y 

The solution for <p{x,t,y) is analogous; in view of this and of the close 
relationship between ^ and $, only \p shall be treated from now on. 

The function \f/(x,0,y) represents initial conditions, and is considered 
as given. To find if/(x,t,0), we use the integral condition (19), and con- 
clude that 

+(x,t,0) = [ +(P v ,t - y,Q)P v a(!/) dy 

- f (P y - l)P y N q N (t - y,0)a(y) dy (22) 

• 

To solve the functional -integral equation (22), we set x = 1 + w, 
and equate coefficients of like powers of w. This yields 

Kit) = f \b„(t - y) + l>n-i(t - y) 

' ' (23) 

e-" yu a(y) dy 4- kM), n ^ 0, 



-C-i)^-^' 



where 



Note that 



Kit) = £(j[)flf(tf>), 

kM = L b " (u)e 1 - A(u) ■ 

iV 

bo(t) = Z?«('.o) - <KU0). 



Let the Laplace transforms of b»( • ), fc„( • ) be b„*(s), k n *(s), respec- 
tively. We obtain a simple recurrence for the &„* by applying the Laplace 
transformation to (23). The recurrence is 

b„* = a n (s) L* + &„_:* - ( ^ J bA + k,*, n ^ 0, (24) 
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where 



A*(s) = f 
Jo 



e- au dA(u), 



a n (s) = A*(s + wy).f 

To find &o*, let x approach 1 in (22) ; then P v (x) goes to 1, and we reach 
the following renewal equation for ba(t)\ 

bo(t) = f bo(t - y)a(u) dy + f / (i y a (t + u) du. (25) 
h .'o 1 — A{u) 

It can be verified that the last term on the right of (25) is just ko(t)', 
upon solving (25) by transforms, we find that b * = /.'o*/[l — A*]. 

It can be seen that b (t) is the density of arrivals at the time t; thus 
b (t) is a familiar function of renewal theory, for which the reader is 
referred to Smith" 1 and the bibliography therein. 

The solution of the recurrence (24) is 



b,* = (x„r i l* - 1 re, n ,) i>»* - -^h 



U X H' 



(2G) 



where 



X = 1, 

1 - a„(s) 



X n = 



a„(s) 



X„_x. 



In particular, the Laplace transform of the density (at t) of arrivals 
finding all trunks busy is given by 



& w * = <i»* = [ e-'qAW) dt = 1 ~ aM v jg aM . (27) 
The generating function of distr jiV(/ — 0)) is 



E\x 



S'il-0) 



fix, t, y) dy 



= f S <hXl - y,0)[P y (x)}" +l ^"[l - A(y)] dy (28) 

+ I * Z ff.(Q* - DIPMY l ~f (y \. dy. 
•' " l — A(y — t) 

t The functions a„(s) are to he distinguished from the constants a„ of Ref. 11, 
which use the same model and notation. The two quantities are related by a„ = 
A*(ny) = a„(0). 
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APPENDIX B 



The Stationary Distribution of [N(t),y(t)} 

We now consider which initial distributions q n (0,u) for (A/(0+), 
,y(04-)} are stationary, i.e., are invariant under the transition probabili- 
ties of the Markov process [N(t - 0),y(t)}, studied in Appendix A. 
Intuitively, since we show in Theorem 3 of Appendix C that a limiting 
distribution exists as t -» » , we expect this limit to give the stationary 
distribution. This is the content of 

Theorem 2: If A(u) has a continuous derivative and mi < co . the 
x,u function 

T,PnPu" +1 - SN "(x) 1 - MU \ u^O, (29) 

n Ml 

generates the unique stationary distribution of \N(0+),y(0)} ; (29) is 
a generating function in x and a probability density in u. The number 
p„ is the equilibrium probability that an arriving customer find n trunks 
busy. 

To show that (29) generates a stationary distribution, it is sufficient 
to prove that the choice of (29) for the initial condition makes each 
q, t (t,0) = p„//ii for all t. This is equivalent to 

Pn 



or to 



g„*(s) = ^, 

sm 



b,*(s) = !± 

S/Ltl 



with 



6„ = E [I )P 



In order to use the recurrence (24) and the formula (27) for q N *, we 
must first calculate the quantities k„* imposed by (29). Now k n (t) is 
the nth binomial moment associated with the generating function 

for yp(x,0,u) given by (29). Thus, k n (t) is associated with 
Mi" 1 E f PnP v J +n - iNn [PM) dA{y). 
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This is equal to 

M1 _1 Z f PnP« +t l+n -* Nn (x) dA(u), 
n Jo 

and so, for n > 0, 

K(t) = mi" 1 Un + />„_, - f '^ J p N \ ^ e-" y "(1A(u). 

The Laplace transform of this is 

i * (iii / N \ \ a n — a,.(s) 

I V"' — V J SMi 



For n = 0, 



W)- 1 -^ 



A-n* = 



Mi 

1 - 4*(«) = 1 - a (g) 

SHi SMi 



We now note that, for these fc n *, the condition b N * = p N /sni implies 
b,* = b„/sni for all lower n. This can be proved by induction from (24). 
We now substitute these /.„* in (27) for q N * ( = b N *). If we divide out 
a factor [1 — o (s)] in the numerator, the first term is l/.s^i ; the general 
term is 



|> + 6 -'-(„-iH 



[1 - ai (s)} •••[!- a n ^(s)][a n - a n (s)) 
(s/ii)ai(s) • • • a n (s) 

Using the recurrence of Ref. 11, 

b n = a„ b n + &„_, - ( J p N \, n > 0, 

we find after much algebra that q N * = p N /sni , which proves the theorem. 
The stationary value ps/n\ for the density q#(t,0) has the following 
physical interpretation: 1/mi is the equilibrium density of arrivals, and 
p N is the chance that such an arrival find all trunks busy. 

The uniqueness of the stationary distribution follows from that of 
the limiting distribution as / —> co. For two distinct stationary distribu- 
tions of necessity give rise to distinct limits, contradicting Theorem 4 
of Appendix C. 

The analog of Theorem 2 for the more general formulation of Appendix 
C is proved by the same form of argument that established Theorem 2, 
with the difference that Laplace-Stieltjes transforms are involved, and 
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that special mention must be made of the "periodic" case, in which 
A(-) is concentrated on a lattice. 

APPENDIX C 

Approach Via Renewal Theory and Regenerative Processes 

This last appendix is a quick sketch of results for general distributions 
A ( ■ ) ; no proofs are given. 

Smith 13 has defined a regenerative stochastic process x(t) as one for 
which there is an event R such that, if R occurs at t, then knowledge of 
x(s) fors<t loses all prognostic value, and the future development of 
x(t) for t > t depends only on the fact that R occurred at t. The points 
at which R occurs are called regeneration points of the process. 

Let R„ denote the event: an arriving customer finds n trunks busy. 
Since the interarrival times form a renewal process, each point in time 
at which R n occurs is a regeneration point of N(t — 0), for all ^ n ^ 
N. In fact, we have already 11 made use of this property of the arrival 
process in const met ing the imbedded Markov chain. 

We are therefore in a position to use Smith's results 13 directly. The 
regenerative property of R» implies that the time intervals between suc- 
cessive occurrences of R n form a renewal process, i.e., a sequence of 
independent, identically distributed variates. To apply the results of 
Ref. 13 we must investigate whether these variates are proper, i.e., 
finite almost surely, and whether they have finite expectations. We 
content ourselves with the following result: 

Theorem 3: Let x m ,„ be the time elapsing from an occurrence of R m 
until the next occurrence of R n . Then 

x m ,„ < co with probability 1, 
and, if the mean interarrival time p x = J x dA < co f then 

E{x m ,„\ < CO. 

We use the following notations: 

iii = the tth interarrival time, i = 1, 2, 3, • • • , 

r,(ra) = the time interval between the (i - l)th and the 
tth occurrences of R,„ , 

k 
U k = ]C u i = tne epoch of the Ath arrival, 
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k 

X k (m) = ^Z Ti(m) = the epoch of the fcth occurrence of the event R m , 

T i = the epoch of the last arrival prior to t and 
after 0. 

The xii all have the common distribution A(u), except for U\ , which 
has G. Also, the r,(m) have a common distribution, except for ri(ra). 

During the interval (T t ,t), the process N(x - 0) forms a pure death 
process whose transition probabilities P m , n (•) are known. Let U(x) be 
the unit step function at and 8 mN the Kronecker delta. The probabil- 
ity that N(t — 0) = n can be represented as 

Pr [N(t - 0) = n) = Elp^^MUiut - t)\ 

+ Z [ Pm+i-s ml/ .n(t - ") du Pr [T, ^ u and N(T t - 0) = ml, 

;i-lSm£AT 

where the measure implicit in the E operation is the joint distribution of 
iV(0+) and u t . With the notations just introduced in mind, it can be 
verified that 

Pr [T t < u and N(T t - 0) = m\ = £ Pr [T, = X k (m) ^ u\ 

k=i 

= f [1 - A(t - u)] d, E Pr [U k £ v and N(U k - 0) = m] 

J k 

= [ [I - A(t - v)] d v J2 Pr \X k (m) ^ v], 

J k 

the series being absolutely convergent. By introducing 

QM) = E Pr [U k ^ I and N(U k - 0) = n\, 

we can write 

Pr [N(t - 0) = n] = E[p Km .i,MU{ui - t)\ 

+ Z / P*+1-W.»0 - W)[l - A (I - K>] dQm(w). 

n— lgm<W ■ 

This representation has been used by Takacs 7 to study lim Pr {N(t — 
0) = m] as / — > =o by methods similar to those used in the proof of 
Theorem 4. 

We can now describe directly some conditions under which Pr \N(t — 
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0) = n] goes to a limiting distribution as t — > 0, independent of initial 
conditions. The result is due essentially to Takacs. 

Theorem 4: If A(u) is not periodic, if ih < <» almost surely, if m = 
E{ui) < oo,i > 1, then 

lim Pr {N(t -0) = n}= S f p n+1 - SmN , n M -pTZ — r~ du - 

This result follows at once from the previous results of this section and 
Theorem 2 of Smith, 13 upon noting that p«,*(w) is a linear combination 
of monotone decreasing functions. From Theorem 3 of Smith 21 there also 

follows 

Theorem 5: If A(u) has period p, if U\ < » almost surely, if mi = 
E[ui\ < oo, i > 1, andO ^ y < p, then 

lim Pr {iV(A:p + y - 0) = n] = 

k-»<*> 

2^ 2. P»+l-»»Y.n(fcp + 2/) p|^ j • 

We now derive and solve equations for the quantities 

Q«(m) = £Pr{tf* ^ w and JV(Z7* - 0) = wi}, 

k 

which occur in the representation for the probability Pr [N(t — 0) = w} . 
Using the generating variable x and the abbreviation 

P tf (aO = 1 + (x - 1)<T™, 
we find that 

2 »" Pr (tft+i ^ < and iV( C/"*+i - 0) - n\ = 



I 

71 = 



2 f f P„ m+1 " imV (.r) dA (i/) d„ Pr { tf* ^ u and tf ( Uu - 0) = to}. 

m JO JO 

The series formed by adding all these equations up on the index k are 
absolutely convergent; hence no additional generating functions are 
needed, and we reach the equation: 

2 x n ZPr [U k ^ t and N(U k - 0) = n} = 

n-0 fc=l 

£ rc" Pr fa ^ t and N(ui - 0) = w} 

»-o (30) 



* ft— u 



+ EJU 

M JO JO 



P y m+1 - im »(x) dA(y) d u 2 Pr { 17* ^ « and 

tf(17* - 0) = to}. 



TRANSITION PROBABILITIES FOR TELEPHONE TRAFFIC 1317 

This is an integral-functional equation for the function 

*0M) = Z> B EPr{tf*g t and N(U k - 0) = n], 

n k 

which is closely related to the function \f/(x,t,0) treated in Appendix A. 
In fact, when ¥ is absolutely continuous in t, then \J/ is its derivative, 
and (22) is similar to (30) in the special case where the density \f/ exists. 
Equation (30) may be solved by the same method as (22), except 
that Laplace-Stieltjes transforms replace the ordinary Laplace integrals 
used for (22). We introduce the following notations: 

Q n (t) for 2 Pi' \Uk £ t and N(U k - 0) = n), 

h 

Bn(t) for f; (fy Qj(t), 

K n {t) for Z ( j ) Pr [ut£t and Nfa - 0) = ;}. 

When each of Q„ , B„ and 7C„ is absolutely continuous, the corresponding 
(lower case) quantities q„(t,0), b ri (t) and k„(l) are the respective deriva- 
tives. Let the respective Laplace-Stieltjes transforms of Q„ , B„ and K„ 
be Q,*, B„* and K„*. Then (30) leads to the recurrence 



B n * = a„(s) Ib,* + £,._!* - ( ^ J B„*j + K n *. 



(31) 



The rest of the solution is in complete analogy with the solution for the 
b„*, Qn* in Appendix A. We find 

BM = Q N (t), 

.V 

B (t) = EQn(0 = *(1,<). 
o 

The function ¥(1,0 satisfies the renewal equation 

¥(1,0 = ( ¥(1,< - y)rfA(y) + 0(1), 
where (? = distr {?/i}. The Laplace-Stieltjes transform of ¥(1,0 is 
,. f e~ 8 ' d(7(0 ^ * 
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The solution of the recurrence (31) is 
B n * = 



ft^Jft'-tKAW-zml 

V 1 - a„(s) \ wLVj- V a M- 



V a.-(«) , 



where the first term of the products is taken to be one. The Q n * are given 
in terms of the B n * by the equation 

Qn* = Z(-ir( n t j ) B ^*- 

j=o \ n / 

In particular, the analog of (27) is 

Bn * = Q N * = P e~ a ' dt Z Pr [U k ^ t and N(U k - 0) = N\ 

JO k 

= [1 - A*(i)r 1 (32) 

oi(s) qi(g) • • • <w) 



/JV\ 1 - ^(s) , (N\ [1 - aM • • • [1 - a w (s)] 



1+ Vl/ axis) + + W ai(«) ■ • • o*(«) 



From the representation of Pr {JV(£ — 0) = n) it can be seen that 
the generating function of N(t — 0) is 

Elx"™) =S{P t A,(0+) (.r)C/(w 1 -0} 

+ E f P.-u 1+ "- anA, (.r)[l - A(f - «)] dQ.(«), 

7. JO 

with Pj(x) = 1 + (•>-■ — l)e -T \ and U the unit step at zero. It. follows 
that the Laplace transform (with respect to of the generating function 
of N(t - 0) is 

P e-"E{x Nit - 0) \ dt = f e~' t E[P l Nm -\x)U{u l - t)) dt 

Jo J o 

+ £ Vd) P ^"P ;/ n+, - i '" v (^)[l - A(y)] dy. 

n JO 

When lim A'{.e A ' ( '" 0) l exists as / — » °o , we can use Abel's theorem 
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for the Laplace transform to evaluate the limits lim Pr {N(t — 0) = n) 
explicitly. As s — » 0, 

lim 8Q n *(s) = lim sG "* (s) 



1 - F n , n *{s) 

1 



■Cf \ X n ,,, j 

But from (32) we find 

Pn 



lim sQ B (s) - ^' 

Ml 



P« 



x <L1 (x) 

Jo 

where p n is the equilibrium probability that an arriving customer finds 
n trunks busy. (These probabilities have been studied in Takacs 7 and 
Benefi," inter alia.) Hence 

E[x n ,n} = mean recurrence time of R„ 

[ xdA(x) 
_ Jo _ Mi 

Pn Pn' 

and, from Theorem 3, 

[l-A(u)) 



mi 



du 



lim Pr {N(t - 0) = n] = Z p m f Pm+is m v „(w) 

lim Eix N( '- (,) \ = Z p„ [ X P v n+l - s »«(z) [1 ~ Mu)] dy, 

with 

P y (x) = 1 + (x - l)e- yu . 
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